Neutrino Bursts from Fanaroff-Riley I Radio Galaxies 
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On the basis of existing observations (at the 4.5 a level) of TeV 7-ray outbursts from the Fanaroff- 
Riley I (FRI) radio galaxy Centaurus A, we estimate the accompanying neutrino flux in a scenario 
where both photons and neutrinos emerge from pion decay. We find a neutrino flux on Earth 
dF„/dE v = 4.5 x 10" 11 (£„/TeV)" 2 TeV" 1 cm" 2 s" 1 , equally spread in flavor as a result of maximal 
mixing. Such a flux will trigger at the IceCube facility about 10 showers/burst, with negligible 
background from atmospheric muons, and primary neutrino energies in excess of 100 TeV. The only 
other FRI radio galaxy observed in the TeV photon energy range at the 4er level is M87. The burst 
nature of this activity is not established; however, we show that the intrinsic neutrino luminosity 
during the active period is the same as the Centaurus A burst. On the assumption that Centaurus 
A typifies the FRI population, we show that IceCube should collect 10 showers (all neutrino flavors) 
in 3 years, attaining a 95% CL sensitivity to the diffuse neutrino flux from FRI radio galaxies in 
one year of observation. 

PACS numbers: 13.85.Tp, 95.85.Ry, 95.30.Cq 



Conventional astronomy spans about 18 decades 
in photon wavelengths, from 10 4 cm radio- waves to 
10" 14 cm 7-rays of GeV energy. Because the universe 
is opaque to photons of TeV energy and above, present 
studies focus on hadrons, neutrinos, and gravitational 
waves as messengers probing the high energy universe. 
The best candidates to serve as messengers in a new as- 
tronomy of the high energy behavior of distant sources 
are neutral particles. This is because the orbit of a 
charged cosmic ray can be substantially bent by the am- 
bient magnetic field of our own galaxy, destroying the 
possibility of locating the source. The most promising 
messenger is the neutrino: it can be copiously produced 
in cosmic beam dumps and can traverse unscathed dense 
astrophysical environments. In this context, we recently 
carried out a study of directional neutrino signals from 
sources in the Galactic plane 0. In this letter we ex- 
tend this project to the extragalactic domain, and de- 
lineate the prospects to identify Fanaroff-Riley I (FRI) 
radio galaxies as high energy neutrino emitters. In our 
discussion it will be seen explicitly how the neutrinos pro- 
vide a signal precisely where photon signals are blocked 
by interactions with the universal radiation backgrounds 
permeating the universe. 

At present, a handful of sources have been established 
as TeV 7-ray emitters ;2]. All of them are nearby BL 
Lac objects characterized by strong rapid variability and 
apparent superluminal motion. These extreme features 
are generally interpreted as a consequence of dissipative 
effects (non-thermal emission) from a relativistic jet ori- 
ented at small angle with respect to the line-of-sight Q. 

There are two principal mechanisms for TeV gamma 
ray production: (i) Electrons undergo bremsstrahlung 
in the magnetic field and/or inverse Compton scattering 
in the ambient photon sea or (ii) the gamma rays are 
directly traced to 7r° decay. Only the second scenario 



can accommodate baryonic cosmic ray production. Since 
such cosmic rays are observed, it is reasonable to assume 
that at least some gamma ray sources operate according 
to the second mechanism. 

In the context of unification models 0, BL Lac objects 
are intrinsically the same as FRI. Though the jet emis- 
sion from FRIs is not strongly Doppler boosted towards 
us (and may even be de-boosted) , in some cases the lack 
of relativistic boosting can be partially compensated by 
proximity to Earth. Therefore, some FRI radio galaxies 
could be detectable in the TeV range. In what follows, 
we examine the consequences of assuming that gamma 
ray emission from FRI originates in 7r° decay, and is nec- 
essarily accompanied by a flux of high energy neutrinos 
emerging from the ir^ population. 

At a distance of 3.4 Mpc, Centaurus A (Cen A) is the 
prototype of FRI galaxies |5| ■ It is the only GeVgamma 
ray source with a confirmed large-inclination jet [6j . Data 
collected in the early 70's with the optical intensity in- 
terferometer operated by Sydney University at Narrabri, 
show a ~ 4.5cr (time average) excess of 7-ray events from 
the direction of Cen A Q • The reported cumulative flux 
averaged over 3 yr of observations is 

F 1 {E 1 > 300 GeV) = 4.4 ± 1.0 x 10" 11 cm" 2 s" 1 . (1) 

Based on a power- law spectrum oc E~ 2 , this corresponds 
to a luminosity at the source of C w 7.5 x 10 40 erg s" 1 . 
It is important to note that since the detector beam did 
not include the radio lobes, it is likely that the source re- 
sides in the compact nucleus of Cen A. The data reveals 
two different periods of heightened activity in the TeV 
range, each of about 1 year duration. Though data at 
the upper end of the spectrum await confirmation, the 
7-ray bandwidth of Cen A could reach energies as high 
as hundreds of TeV. There are specific hints in this di- 
rection in data reported during the 1980's and 1990's by 
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the Buckland [8j and JANZOS 9] collaborations. Several 
items can be noted with respect to these measurements: 
(i) The data are consistent with bursts of approximately 
one year duration every decade, (ii) Both experiments 
reveal an upper cutoff of ~ 200 TeV, which can be as- 
cribed to absorption on the background radiation fields. 
(in) The reported flux in both cases is about 2 orders 
of magnitude larger than that measured at the Narrabri 
Observatory; the chance probabilities are 6 x 10~ 3 and 
2 x 10 -2 for Buckland and JANZOS, respectively. Al- 
though these data are not at the discovery level, the ob- 
servations of 100- TeV photons at the two facilities corrob- 
orate one another, (iv) The emission of 100-TeV gamma 
rays is supportive of their origin in pion production and 
decay processes, as opposed to synchrotron and inverse 
Compton scattering jlOj . 

Additional relevant information about the source is 
contained in these observations. The ~ 1-year duration of 
the observed burst at Narrabri implies a coherent region 
of activity r ~ 0.3 pc. Therefore, the recent discovery 
of a sub-parsec radio counterjet in the nucleus of Cen A 
is of interest As a conservative estimate we adopt 

the flux in Eq. JIJ as the photon excess characterizing the 
Cen A burst, which corresponds to an integrated burst 
energy E hulst = 2.3 x 10 48 erg. 

A high-energy neutrino flux emerges if the charged 
pion decay length is smaller than the pion interaction 
length in the source region. The latter, as well as the 
pion production rate, depends on the gaseous hadronic 
particle density. An estimate of the target nucleons can 
be obtained by assuming that the gas density n near the 
accretion disk of the supermassive black hole engine of 
Cen A is similar to that near the center of our Galaxy. 
This in turn is arrived at by assuming that the gas/star 
density ratio near the center is the same as the one in 
the Galactic disk, yielding n ss 10 6 cm -3 . Our estimate 
is most likely a lower bound on the gas density, since 
the Cen A black hole mass Af B H ~ 2 x 10 8 111 i s 
about 2 orders of magnitude greater than Sgr A* [lj, [l4( . 
With this in mind, the corresponding pion mean free path 
(ncr-n-jv) -1 w 6.5 pc is much greater than the charged 
pion decay length, which at 100 TeV is 5.6 x 10 8 cm. 
(Here avN ~ 50 mb is the pion-nucleon cross section.) 
On the other hand, the mean free path for collision of 
the ultrarelativistic accelerated protons on the gas is 
{rujpN)^ 1 sa 4.6 pc, where <t p n ~ 70 mb [l^. Assum- 
ing no significant deflections on the magnetic field, this 
implies a probability of interaction in the coherent region 
p = r n (T p N ~ 7%. The infalling mass Mj n f a n required to 
power the burst is found through the relation 



sufficient pion production. 

Since 7r n 's, 7r + 's, and 7r _, s are made in equal num- 
bers, one expects two photons, two ^ e 's, and four v^s 
per 7T°. On average, the photons carry one- half of the 
energy of the pion, and the neutrinos carry one-quarter. 
The energy-bins dE scale with these fractions, and we 
arrive at 
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for the fluxes at the source, where it denotes any one of 
the three pion charge-states. Terrestrial experiments (see 
e.g. (17) 1 have shown that and v T are maximally mixed 
with a mass-squared difference ~ 10 _3 eV 2 . This together 
with the known smallness of |(^ e |f3)| 2 , implies that the 
fu's will partition themselves equally between v^'s and 
u t 'b on lengths large compared to the oscillation length 
A osc ~ 1.5 x 10~ 3 (E u /PeV) pc. Here v z ~ (u„ + u T )/Vl 
is the third neutrino eigenstate. From these remarks, one 
finds a nearly identical flux for each of the three neutrino 
flavors (j = e, fj,, r), which is equal to 
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On the assumption of an E~ 2 spectrum, we fix the nor- 
malization using the cumulative number flux of Eq. 
with an upper cutoff at 3 TeV Q ■ From Eq. (0J we then 
obtain the neutrino flux on Earth 

The observation of 100- TeV photons supports the extrap- 
olation of the spectrum in Eq. JSJ to those high energies. 

The most promising prospect for detection of such a 
low flux is the IceCube facility under construction at the 
South Pole 18]. For energies above - 100 TeV, IceCube 
will provide 4-7r detection of neutrinos since the atmo- 
spheric muon flux is negligible. Moreover, at these en- 
ergies the telescope will be able to resolve directionality 
in all neutrino flavors. The total number of showers in a 
single burst of duration t ~ 1 yr is then 
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where G is Newton's constant and r s is the Schwarzschild 
radius. This gives M infaU » 3.6 x 10" 5 M . This is 
about 100 times larger than the estimated infall rate for 
the Galactic black hole 0] . Because of the burst nature, 
coupled with the much larger mass for the Cen A black 
hole, we take this as a reasonable accretion rate, allowing 



where Na is Avogadro's number, V c s » 2 km 3 is the ef- 
fective volume of ice with density p, and = 6.78 x 
10~ 35 (E u /TeV) - 363 cm 2 is the charged current neutrino- 
nucleon cross section . By substituting the flux given 
in Eq.©, we find about 10 showers per burst, with pri- 
mary neutrino energy > 100 TeV. Thus we arrive at the 
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first prediction of this paper: during a burst similar to 
the 4-5cr event detected by the Narrabri Observatory, Ice- 
Cube will detect about one 100-TeV neutrino per month 
pointing to Gen A. This signal persists even in the ab- 
sence of the emitted photons which are absorbed on the 
infrared and microwave radiation backgrounds. 

Located in the Northern hemisphere at a distance of 
16 Mpc, M87 is the FRI with the brightest optical jet |2C| . 
It shows most of the characteristics of BL Lac objects, 
with the jet oriented at 30° — 35° to our line of sight. 
Data taken during 1998-1999 with the HEGRA stereo- 
scopic system of 5 imaging atmospheric Cerenkov tele- 
scopes show an excess of photons from the direction of 
M87 with significance level 4.1a above background |2lj . 
The observed cumulative flux is 



F 1 {E 1 > 730 GeV) 



0.96±0.23xl(T 12 cm" 2 s" 1 



(7) 



The data can be fit with a power law dF 1 jdE 1 
with a — 2.9± 0.8 s tat ± 0.08 sys t. The large uncertainty in 
the spectral index is not, however, reflected in obtaining 
the source luminosity: with an upper cutoff of 5 TeV, 
there is about a 25% variation as a varies between 2 and 
2.9. Hereafter, we take a = 2 to reflect the conventional 
Fermi engine emission spectrum. This corresponds to a 
7-ray luminosity at the source £ w 6.8 x 10 40 erg s _1 . 
This is remarkably close to the 7-ray source luminosity 
of the Cen A burst obtained above. 

By duplicating for M87 our Cen A discussion, we ob- 
tain 
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A rough estimate of the event rate at IceCube can 
be obtained following the analysis of our previous pa- 
per pj. Using the terrestrial transmission probability 
given in [23, and a conservative effective area of 1 km 2 , 
we find about 2 events/yr during a period of TeV gamma 
ray activity. This is similar to the atmospheric back- 
ground. A more definitive assessment of sensitivity to 
such a signal will await further refinement of angular and 
energy resolutions via improved knowledge of the detec- 
tor response. Additionally, the temporal profile of future 
gamma ray observation will allow a better delineation of 
the signal. 

Although there are no other nearby FRIs of this magni- 
tude which can potentially be detected as point sources, 
one can integrate over the estimated FRI population out 
to the horizon to obtain a prediction for the diffuse neu- 
trino flux. This quantity is given by 
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where 1Z ~ 1 horizon ~ 3 Gpc, n FRI ~8x 10 4 Gpc~ 3 is 
the number density of FRI [23, and £„ = (dNo/dE^/r 
is an average neutrino luminosity (all flavors) of FRI ra- 
dio galaxies. Here dN /dE u is the differential injection 
spectrum for a single burst and r ~ 10 yr is the period 



between outbursts. Although the burst nature of the 
HEGRA observation is not established, in what follows 
we conservatively assume that the photon excess from 
M87 is localized to an outburst period. Since both Cen 
A and M87 have almost the same luminosity, we adopt 
this as the average for FRI radio galaxies. Thus, the 
quantity dNo/dE^ can be obtained from Eq. (J5J), 
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where Nf = 3 is the number of neutrino flavors and d 
is the distance to the source. For long time averages we 
find 
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The diffuse flux observed on Earth then follows from 
Eq. © 



El J„ ~ 1.2 x 10" 11 TeV cm~ 2 s" 1 sr" 1 



(12) 



Note that this flux is about a factor of 2 smaller than 
the Waxman-Bahcall [24| upper limit on the intensity 
of neutrinos produced in sources which also emit bary- 
onic cosmic rays. After 3 years of observation, the 
90% CL sensitivity (corresponding to 2.44 events) of Ice- 
Cube to the diffuse muon neutrino flux is E 2 J„ ~ 
3 x 10~ 12 TeV cm" 2 s _1 sr" 1 0. We arrive then at 
the second prediction of this work: If the observations of 
Cen A at Narrabri and HEGRA measurements of M87 
characterize the outburst of FRI radio galaxies, then Ice- 
Cube will collect about 10 neutrino events (all flavors) in 
three years. 

It seems worthwhile to briefly examine the implications 
of neutrino bursts of magnitude corresponding to the 7- 
ray intensities reported by the Buckland and JANZOS 
collaborations, with full knowledge that these observa- 
tions arc not at the discovery level. Such a powerful 
burst implies a neutrino flux on Earth about 2 orders 
of magnitude greater than that in Eq. JSJ. In addition 
to vastly increasing the event rate at IceCube, the 
flux should be detected by 0.1 km 2 neutrino telescopes 
under construction (ANTARES [25|, NESTOR [H) and 
a planned 1 km 3 facility (2?|], all in the Mediterranean. 
The location of ANTARES at 43° North provides a 5cr 
discovery sensitivity of 1.2 x lO" 10 {E u /TeV)- 2 TeV" 1 
cm -2 s _1 (for muon neutrinos) in the direction of Cen 
A [28j. We find that in the course of ~ 10 years of obser- 
vation, ANTARES and NESTOR will measure a neutrino 
burst an order of magnitude larger than the 5<r detector 
sensitivity. 

In summary, we have analyzed the possibility of detect- 
ing the neutrino counterparts of various TeV 7-ray obser- 
vations in a model where gamma rays originate through 
7r°-decay at the source. We have found that Icecube will 
attain sensitivity to observe neutrino bursts from Cen A. 
Moreover, if the 7-ray observations with the Narrabri Ob- 
servatory and the HEGRA Cerenkov telescopes charac- 
terize the emission behavior of the FRI population, then 



4 



IceCube should observe the diffuse flux with statis- 
tical significance in several years. In light of this, we 
conclude that observations at future neutrino telescopes 
will permit a major advance in discriminating between 
high energy astrophysical processes. 
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